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ABSTRACT 


A  sinusoidal  surface  profile  is  used  to  study  by  an  exact  method 
the  effect  of  water  waves  on  an  electromagnetic  field  propagating  downwards 
from  the  surface.  It  is  assumed  that  the  magnetic  field  is  directed  parallel 
to  the  surface  corrugations.  The  results  are  presented  graphically. 

Comparisons  are  made  between  these  results  and  those  obtained 
using  an  approximate  method  of  Wait. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 


Water-Wave  Effects  on  Radio  Wave  Propagation  in  the  Ocean 


I.  METHOD  OF  SOLUTION 

The  problem  under  consideration  is  the  evaluation  of  the  electro¬ 
magnetic  field  components  under  the  ocean  when  the  field  sufficiently  far 
above  the  ocean  is  essentially  uniform  and  horizontal.  In  a  Cartesian 
coordinate  system  with  its  z  axis  pointing  vertically,  the  ocean  surface 
is  assumed  to  be  given  by  the  equation 

z+acosKx  =  0.  (1) 

where  k  =  2tt/l  and  L  is  the  water  wavelength.  The  magnetic  field  H 
far  above  the  surface  is  assumed  to  lie  in  the  y-direction,  and  so  the 
symmetry  of  the  problem  forces  it  to  remain  so  everywhere.  Since,  in 
addition,  neither  the  far  field  nor  the  surface  shape  is  a  function  of  y,  then 
no  field  quantity  will  depend  on  y. 

Now  the  electromagnetic  wavelengths  of  practical  interest  are  very 
long  compared  with  the  water  wavelength.  Therefore  the  quasi- static 
assumption  that  displacement  currents  are  negligible  can  safely  be  made. 
Then  in  the  (non-conducting)  air,  H  satisfies  the  equation  V  x  H  =  0, 

V.  H  =  0.  But  since  H  =  yH^(x,z),  then  the  only  possible  solution  in  the 
air  is  that  H  is  a  constant  evervwhere. 

y 

2  2 

In  the  ocean,  H  satisfies  the  quasi- static  equation  (V  +  k  )H  =  0, 

2  ^ 

where  k  =  itDOp^.  (Time  dependence  is  assumed  to  be  as  exp|-rJUt|,  O 
is  the  ocean  conductivity  and  |Ji^  is  the  permeability  of  free  space).  Thus 
since  the  solution  must  also  be  periodic  in  x,  the  general  solution,  if  it 
converges,  for  H  is 


1 


(2) 


H  (x,  z)  =  7,  b  exp(-ik  z)cosnKx. 
y  ^  n  -n 

^  n=o 


n 


(The  symmetry  of  the  problem  about  the  plane  x  =  0  excludes  the 
possibility  of  additional  terms  in  sinnKx).  Here  k  is  defined  as 


n 


n 


/  2  2  r 

=  vk  -  n  K 


.  /  2  2  A 

=  lyn  K  -  k 


(3) 


and  the  negative  sign  in  the  exponent  in  (2)  ensures  that  each  mode  decays 
exponentially  as  z  goes  negative. 

The  corresponding  expression  for  the  electric  field  E  (x,  z)  is  given 
by  E  =  (VxH)/a  =  (-xB  H^/Bz  +  z  B  H^/Bx)/a  ,  or 

00 

E  =  —  /.  b  Ixik  cosnKx  -  z  nK  sin  n  K  x  lexp( -ik  z)  (4) 
a  ^  n*  n  f  n  '  '  ' 

n=o 


Solving  the  problem  involves  finding  the  values  of  the  unknown  b 


n 


coefficients  in  (2)  and  (4)  by  means  of  the  boundary  condition  that  H  is 
continuous  at  the  surface.  Thus,  assuming  the  uniform  field  above  the 
surface  is  normalized  to  unity,  one  can  formally  write,  using  (1)  and  (2), 


^  ^n  cosKx[cosnKx  =  1. 


(5) 


n=o 


By  truncating  the  series  on  the  left  to  N  terms  and  enforcing  the 
equality  at  N  points  in  the  range  0  ^  x  ^  E  one  can  obtain  N  equations  in 
the  first  N  unknown  b^  coefficients.  These  equations  can  then  be  solved 
by  the  usual  methods.  Unfortunately,  to  obtain  (5),  one  must  make  the 
Rayleigh  assumption,  which  is  that  the  downward  -  going  wave  expansion  (2) 


2 


is  valid  not  only  for  z  <  -  a  but  also  in  the  strip  -a  ^  z  ^  a.  It  can  be  shown 
[  1,  Z  ]  that  if  the  normalized  wave  height  a/L  is  greater  than  0.  713  then 
(2)  is  invalid  there  and  (5)  diverges.  Thus  the  b^  found  by  truncating  the 
series  in  (5)  to  N  terms  will  not  in  general  converge  to  their  correct  values 
as  N  is  increased. 

However,  by  interpreting  (5)  in  a  more  general  sense,  one  can  still 
use  it  as  a  basis  for  finding  the  b^  even  though  the  series  diverges  [  3,  4  ]  . 
That  is,  one  regards  the  left  side  of  (5)  as  a  generalized  function  which 
can  be  equated  to  the  actual  field  on  the  surface  only  indirectly  via  a  complete 
set  of  sufficiently  smooth  test  functions.  In  this  sense  the  downward  going 
wave  expansion  will  remain  valid  on  the  surface  for  a  much  larger  range  of 
a  than  in  the  conventional  point-by-point  sense. 

Thus  although  (5)  may  be  incorrect  as  it  stands,  the  result  of  multi¬ 
plying  each  side  by  (K/2n)cosmKx  and  integrating  with  respect  to  x  over 
the  range  0  to  2tt/k  is  the  equation 
00 

y  A  b  =  6  (6) 

^  mn  n  mo 
n=o 


where  6  =  1  if  m  =  0  and  is  zero  otherwise,  and  A  is 

mo  mn 


Ztt/k 

K  p  ' 

A  =  exp(i  k  acosKx)  cosnKx  cosmKx  dx 


^|in+mj  +  (k  a)[/2 

»  n+m  n  n-m  n  > 


(7) 


Here  the  chosen  te  st  functions  are  the  set  (k/Ztt)  cos  m  Kx,  (k/Zh)  sin  m  k  x, 
of  which  only  the  former  are  necessary  for  representing  the  even  function 
of  X  which  is  H  ,  and  the  series  in  (6)  converges  over  a  much  larger  range 


3 


of  a  than  does  the  series  in  (5).  In  (7),  J  ,  (k  a)  is  the  Bessel  function 
of  the  first  kind  [5]  of  order  n+m  and  argument  k^a. 

The  evaluation  of  the  is  now  carried  out  in  a  straightforward  way 
by  truncating  the  series  and  solving  the  resulting  finite  set  of  equations. 

To  obtain  the  actual  field  quantities  in  the  ocean  from  the  b^  is  a 
simple  matter  in  the  region  z<  -  a  where  the  downward  going  wave  expan¬ 
sion  is  known  to  converge.  One  simply  substitutes  the  numerical  values  of 

the  b  into  (2)  for  H  and  into  (4)  for  E.  However,  in  the  region  -a^z^a, 
n  y  - 

the  series  in  (2)  and  (4)  will  not  in  general  converge,  and  so  cannot  be  used 
directly. 

A  general  method  which  is  applicable  for  evaluating  the  field  quantities 
anywhere  beneath  the  surface  is  to  take  that  function  of  x  which  is  the  field 
quantity  along  the  line  z  =  ^  -acosKx  -  d  and  expand  it  in  the  test 

function  set.  In  this  case,  for  H  for  example,  an  even  function  of  x,  the 

y 

expansion  is 


H  (x, 

y 


00 

n=o 


(8) 


But  the  generalized  function  for  this  same  H^(x,  z^)  is  given  by  (2)  with  z 
replaced  by  z^(x).  Thus  the  c^(d)  can  be  evaluated  by  multiplying  the 
right  sides  of  both  (2)  and  (8)  by  (k/2tt)  cos  m  K  x  and  integrating  from  0 
to  2tt/k.  The  result  is 

00 

c  (d)  =  e  y  A  exp(ik  d)b  ,  (9) 

n  n  ^  nm  m  m 

m  =o 


where  e  =  1  if  n  =  0  and  is  2  otherwise,  and  both  the  b  and  the  A 

n  m  nm 

[defined  by  (7)  ]  are  already  known. 
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Thus  by  using  (8)  and  (9)  one  evaluates  H  on  contours  which  are  the 
same  shape  as  the  surface  contour  but  at  an  arbitrary  depth  d  below  it. 
(When  d  =  0  the  contour  lies  at  the  surface,  on  which  H  =1.  Thus  from 

y 

(8),  the  c^(0)  are  given  by  c^(0)  =  6^^  and  so  (9)  reduces  to  (6),  the 
equation  determining  the  b^.  ) 

Similarly,  by  expressing  E  on  the  shifted  surface- shaped  contour  as 
00  00 

E  =x  /,  a  (d)cosnKx-z/,  s  (d)sinnKx 
^  z-/  n  n'  ' 

n=o  n=l 


one  finds 


a  (d) 
n 


00 

e  y.  A  exp(ik  d)  ik  b  /a 
n^  nm  m  mm 

m  =o 


and 


s  (d) 
n 


00 


m  =  l 


B 

nm 


exp(ik  d) 
m 


m  K  b  /  O’ 

m' 


where  A  is  defined  bv  (7)  and 
nm  ^ 


B 

mn 


.  .n-m  _  ..  ^ .  .n+m  _ 

n-m  n  n+m 


(k^a)]/2. 


It  should  be  noted  that  although  the  method  described  above  can  be 
expected  to  converge  numerically  to  the  exact  values  for  the  coefficients 
bm>  it  will  not  in  general  be  the  case  that  |b  |  goes  to  zero  as  m  grows 

indefinitely,  for  the  series  in  (2)  does  not  in  general  converge  when  z  =  0. 
Since  (2)  does  converge  for  z  ^ -a,  however,  one  should  find  that  the  modi¬ 
fied  coefficients  b  exp(ik  a)  do  go  to  zero  in  absolute  value  as  m  grows 

m  m 

indefinitely. 
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II.  RESULTS 


Approximations  of  the  coefficients  were  obtained  by  inverting 

the  following  truncated  version  of  (6) 


N 

y  A  b  =6  (m  =  0,  1,  .  .  .  ,  N) 
^  mn  n  mo 
n=o 


(10) 


for  various  values  of  N  and  of  the  problem  parameters  6/l  and  a/L. 

(K  =  Ztt/l  and  k  =  (1  +  i)/6,  where  6  =  Jz/^)  is  the  skin  depth  in  the 

water.  )  As  a  check  on  the  accuracy  of  the  inversion,  the  b  values  found 

^  n 

from  it  were  substituted  back  in  the  left  side  of  (10)  and  the  result  compared 

with  6  .  As  a  check  on  the  complete  approximation,  the  sum 

mo 


S  =  Yh  exp(ik  a) 

Lj  n  ^  n 

was  computed  using  the  approximate  b^  and  compared  with  the  correct 
value,  given  by  (2)  with  z  =  -a  and  x  =  0,  of  unity.  All  computations  were 
performed  on  the  IBM  360  computer  using  single  precision. 

It  was  found  that  the  method  converged,  as  N  increased,  to  give  stable 
values  for  the  b  when  a/L  was  about  0.  3  or  less.  The  rapidity  of  the 

m  I-  / 

convergence  was  greatest  for  the  smaller  values  of  a/L.  The  effect  of 
6/L  on  the  convergence  was  less  marked,  but  was  in  the  direction  that 
convergence  was  less  rapid  for  the  smaller  6/L  values.  A  typical  result 
is  that  with  6/L  =  1,  a/L  =  0.  15  and  N  =  9,  the  value  of  S- 1  was  found 

_  3 

to  be,  in  absolute  value,  less  than  10 

When  values  of  a/L  larger  than  0.  3  were  used,  the  convergence  was 

slow  and  the  b  did  not  settle  to  stable  values.  However,  the  inversion 
m 

check  showed  that  this  behavior  was  always  accompanied  by  poor  inversion 
accuracy.  This  indicates  that  the  set  of  equations  (10)  becomes  poorly 
conditioned  when  a/L  exceeds  0.  3,  and  then  single  precision  is  insufficient 
to  solve  them  accurately. 
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Some  specific  results  for  the  underwater  field  components  are  given 
in  Figs.  1  and  2.  In  Fig.  1,  the  horizontal  field  components  and 
in  the  plane  z  =  -a  are  plotted  as  a  function  of  x  in  magnitude  and  phase 
over  a  half-period  (O^x^L/Z)  for  various  values  of  6/l.  The  wave  height 
parameter  a/L  is  constant  at  0.  15  (giving  a  maximum  wave  slope  of  about 
1)  for  all  curves,  and  the  field  quantities  are  normalized  with  respect  to 
the  fields  that  would  exist  at  depth  a  beneath  a  plane  surface. 


The  points  plotted  in  Fig.  1  are  the  result  of  using  Wait*s  [6]  approx¬ 
imation  to  calculate  the  same  field  quantities.  This  method  assumes  that 
for  sufficiently  shallow  and  long  water  waves,  the  field  propagates  down  to 
the  level  z  =  -a  essentially  as  a  plane  wave.  Thus  it  would  be  expected  to 
give  good  results  for  small  6/L  and  small  a/L.  The  extent  to  which  the 
approximation  deteriorates  when  a/L  is  0.  15  is  indicated  by  the  closeness 
of  the  points  to  the  continuous  lines.  (In  Wait*s  approximation,  the  nor¬ 
malized  quantities  H  and  E  are  identical,  so  that  only  one  curve  of 
^  y  X 

points  exists  for  amplitude  and  one  for  phase.  ) 


At  a  sufficiently  great  depth,  only  the  first  terms  in  the  downward¬ 
going  wave  expansions  (2)  and  (4)  remain  significant.  Under  this  condition, 
the  ratio  of  the  actual  field  to  the  field  that  would  exist  at  the  same  depth 

beneath  a  plane  surface  is  simply  b  /l  for  both  H  and  E  ,  since  b  =1 

^  o  y  X  o 

when  a/L  =  0.  It  was  found  that  |b^|  is  always  greater  than  or  equal  to 
unity,  and  so  the  amplitude  departure  can  be  represented  unambiguously 
by  the  quantity  |b^  -  1  ].  The  phase  departure  is  just  the  phase  of  b^.  These 
quantities  are  plotted  in  Fig.  2  as  a  function  of  a/L  for  various  values  of 

6/l. 


Since  b^  also  equals  the  ratio  of  the  average  horizontal  field  to  the 
field  at  the  same  depth  below  a  plane  surface,  Fig.  2  also  shows  the  way 
the  average  horizontal  field  components  depend  on  a/L  and  6/L  for  any 
z  ^  a. 
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III.  CONCLUSIONS 


The  mathematical  technique  described  here  (a  generalized- 
function  interpretation  of  the  Rayleigh  assumption)  appears  to  be  accurate 
and  to  converge  rapidly  for  values  of  a/L  for  which  the  conventional 
interpretation  of  the  Rayleigh  assumption  is  known  to  be  invalid.  Using 
this  technique,  one  can  evaluate  the  electromagnetic  field  at  any  point 
beneath  the  surface  for  all  water  wave  heights  and  lengths  of  practical 
interest. 


REFERENCES 

1.  Petit,  R.  and  Cadilhac,  M.  ,  "Sur  la  diffraction  d*une  onde  plane  par 
une  re'^seau  infiniment  conducteur",  CR  Acad.  Sci.  ,  1966,  Z6ZB, 
pp.  468-71. 

Z.  Millar,  R.  F.  ,  "On  the  Rayleigh  Assumption  in  Scattering  by  a 

Periodic  Surface",  Proc.  Camb.  Phil.  Soc.  ,  1969,  65,  pp.  773-91. 

3.  Burrows,  M.  L.  ,  "Equivalence  of  the  Rayleigh  Solution  and  the 
Extended-Boundary-Condition  Solution  for  Scattering  Problems", 
Electron.  Lett.,  1969,  5,  pp.  Z77-8. 

4.  Burrows,  M.  L.  ,  "Example  of  the  Generalized- Function  Validity 
of  the  Rayleigh  Hypothesis",  Submitted  to  Electron.  Lett. 

5.  Jahnke,  E.  and  Emde,  F.  ,  Tables  of  Functions  (Dover,  1945),  p.  149. 

6.  Wait,  J.  R.  ,  "The  Calculation  of  the  Field  in  a  Homogeneous  Conductor 
with  a  Wavy  Surface" ,  Proc.  I.  R.  E.  ,  1959,  pp-  1155-6. 


8 


S/L  =  0.1  1  2.0  L  S/L 


aanindww 


9 


Fig.  1.  The  amplitude  and  phase  of  the  horizontal  field  components  and  Hy  in  the  plane 
of  the  wave  troughs  as  a  function  of  x/L  for  a/L  =  0.  15  and  various  values  of  6/L.  The 
fields  are  normalized  with  respect  to  the  fields  at  the  same  depth  (z  =  -a)  beneath  a  plane 
surface.  The  plotted  points  were  calculated  using  Wait's  approximation. 


a/L 


Fig.  2.  The  amplitude  and  phase  departures  of  the  horizontal  field  at 
•’great"  depth  (or  of  the  average  horizontal  field  at  any  depth  z  ^  -a) 
from  those  of  the  field  at  the  same  depth  below  a  plane  surface. 
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